Terahertz absorption by excitonic polaritons 
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We use linear response theory in order to compute the light absorption spectrum, in the terahertz 
band, of a polariton system composed by excitons in a quantum dot very strongly coupled to 
the lowest photon mode of a thin micropillar. In a thermalized (Bose condensed) system at low 
temperatures, the spectral function shows a peak associated to a ls — 2p exciton transition, enhanced 
by polariton effects. On the other hand, in a non-equilibrium system absorption is peaked at low 
energies. Thus, a measurement of terahertz absorption could give an indication of the degree of 
thermalization in the polariton system. 



The strong coupling regime in the interaction between 
a confined photon mode and electron-hole pairs in semi- 
conductor nanodcvices has been demonstrated recently 
[1]. The quasiparticles, so called polaritons [21 13], which 
are roughly half excitons and half photons, offer very in- 
teresting possibilities, such as, for example, a new lasing 
mechanism (polariton lasing) based on their quasibosonic 
nature |4j, with pumping threshold (related to ground- 
state occupation) two orders of magnitude lower than 
ordinary (photon) lasing in the same devices [5j, and op- 
eration at ambient temperatures [6]. 

In the present paper, we focus on the linear response 
of a model polariton system to terahertz radiation. The 
first motivation to carry on such a study is the intuitive 
idea that the interaction with the confined photon mode 
reinforces coherence of the excitonic subsystem and, thus, 
may reinforce the collective response of the excitons to 
the terahertz probe. This may result in a semiconductor 
version of the Giant Dipole Resonances (GDR), a phe- 
nomenon widely studied in nuclei [7] and electron clusters 
[S^ , with the possibility of controlling the position and in- 
tensity of the resonance by varying parameters such as 
the pumping rate or the photon-exciton detuning. 

The second good reason to study terahertz absortion 
by excitonic polaritons is that it has proven to be very 
useful in order to observe exciton formation dynamics 
in quantum wells and bulk systems [10]. We think, 
the available experimental techniques should be able to 
measure the degree of thermalization of the polariton 
system, not only under stationary conditions ^Ij, but 
in the pumped regime as well [12]. Indicators following 
from interband emission alone are not enough because 
the main qualitative features (population of the lowest 
polariton state, behavior of the second order coherence 
function, etc) can be reproduced also from dynamical 
equations, without any thermalization mechanisms, both 
in the pumped 13J and in the stationary regimes [14j . 

Below, we compute terahertz absorption in two ex- 
treme situations. One is a Bose condensed state at low 
temperatures, described in a canonical ensemble. It is as- 
sumed that only the ground state has a significant occu- 



pation probability. The spectral function shows a GDR- 
like peak, whose position grows with the polariton num- 
ber. The second case corresponds to a polariton system 
in a non-equilibrium stationary state (result of a balance 
between pumping and losses), with occupation probabil- 
ities that can not be described by a Gibbs distribution. 
The terahertz spectral function gets a completely differ- 
ent shape with a central peak at near zero energy. In- 
termediate situations would interpolate between the two 
extremes, and a measurement of the response in real sys- 
tems would indicate their degree of thermalization. Let 
us stress, however, that calculations are carried on in a 
model with very strong light-matter coupling, in which 
Rabi splitting overpasses the exciton binding energy. 

(i) Ground-state response of non-interacting 
polaritons 

In order to get a preliminary estimate of the absorption 
spectrum, we first consider the ground-state response of 
non-interacting polaritons. We assume the system is 
in a Bose-condensed state, with Npoi polaritons occu- 
pying a single state. Intraband absorption is described 
by the dipole operator acting only on the exciton func- 
tions. The absorption probability is then proportional 
to \a diol'^Npoi, where a is the Hopfield coefficient [2] 
(that is, amplitude of the exciton ground state func- 
tion in the polariton function), and dio is the intra- 
band dipole matrix element between ground-state exci- 
ton and an excited-state function. The latter is supposed 
to concentrate the oscillator strength for dipole transi- 
tions. Notice that the absorption probability increases 
with the number of polaritons in the ground state. The 
peak position, on the other hand, should be roughly con- 
stant, equal to the energy difference between the exciton 
ground- and excited states. 

Finite, but low, temperatures, should lead to similar 
results. In a grand canonical description, on the other 
hand, which is more natural for the polariton system, 
sectors with polariton number near the mean value will 
contribute also to the spectral function with relatively 
high weights. The effects of polariton-polariton interac- 
tions is considered in the next paragraph. 
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1: (Color online) (a) The ground-state spectral func- 
Eq. fl]), for various Npoi numbers, (b) The non- 



FIG. 

tion, 

equilibrium spectral function, Eq. ( |12[ | , for pumping rates (in 
ps~^) corresponding to mean polariton number in the steady 
state, (Npol), in the interval (1,10). The detuning parameter 
is A = —3 meV. 



(ii) Ground-state response of interacting polari- 
tons 

Polariton-polariton interactions come from residual 
Coulomb interactions between excitons. Instead of us- 
ing a phenomenological approach, we start from a model 
in which Coulomb interactions are treated exactly, and 
the fermionic degrees of freedom are explicit. There is 
a finite number (10) of single-particle states for electron 
and holes, and a single photon mode. Saturation effects 
due to Fermi statistics are seen when the polariton num- 
ber is around (or greater than) 10. A detailed description 
can be found elsewhere [TJ- 

The very-low temperature (ground-state) response of 
the A'poi-polariton system is contained in the spectral 
function: 
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where matrix elements, {I\d\J), of the intraband dipole 



operator, d '-^ 



) (where r^''-' and f^*^' are. 



Pi 
Q 

O 0.005 



Npol = 4 



#- ® A = — 3 meV 
•- • A = + 3 meV 



9 Npol = 5 

/ "~ ^ 

/ Npol = 1 



Npol= 10 



8 10 12 14 16 18 20 

GDR position (meV) 

FIG. 2: (Color online) Intensity (that is, dipole matrix ele- 
ments squared) and position of the GDR peak for two different 
values of the detuning, A. Each dot corresponds to a given 



be computed. | J) is the ground state function of the Npoi- 
polariton system, and \I) are excited states. Fq — 0.1 
meY/h is a phenomenological damping parameter, and 
LOij ~ {El — Ej)/h - the transition frequencies. 

In our model, wave functions are constructed as linear 
combinations; 
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where Se and Sh are Slater determinants for electrons 
and holes, with electron and hole numbers and TV/j, 
respectively, and n is the number of photons in the con- 
fined mode. Functions entering the combination preserve 
the polariton number: 



pol 



N, 



(3) 



and the total (envelope) angular momentum projection 
along the cavity axis (we assume a circular section): 



(4) 



respectively, the hole and electron position vectors) shall 



The ground-state function, |J), has L = 0, whereas |/) 
are L — 1 functions. 

We show in Fig. [l] (a) the spectral function for dif- 
ferent polariton numbers and detuning A = — 3 meV. In 
the model, the parameter A measures the photon energy 
with respect to the nominal band gap, not the photon- 
cxciton detuning. A = — 3 meV approximately corre- 
sponds to resonant conditions. 

The GDRs can be identified as the dominant peaks in 
these curves. The peak position monotonously increases 
with increasing polariton number. This can be under- 
stood on intuitive grounds. The mass of the electron (or 
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hole) cloud is to Npairs, and the Hooke coefRcient for 
the force acting between clouds is fc ~ N^^^^^. Then, the 

excitation energy of the dipole mode is tvjj ~ yjk/m ~ 
\/ Npairs ^ \/ Npol ■ The maximum intensity, on the other 
hand, has a non-trivial dependence on Np^i , a kind of sat- 
uration effect is observed. The intensity first increases, 
as in the non-interacting case, but then, after reaching 
a maximum value, decays. These dependences are illus- 
trated in Fig. [2] where the case A = +3 mcV, corre- 
sponding to an enhanced excitonic component of polari- 
tons, is also shown. In this positive detuning situation, 
the absorption probability rises because the Hopfield pa- 
rameter a increases. 

In spite of the fact that calculations are performed in 
a particular model, we expect that the statement about 
the existence of a peak in the absorption spectrum at 
relatively high excitation energies (of the order of the 
exciton Is — 2p transition), whose intensity increases at 
least for polariton numbers well below saturation values, 
is general enough, and could be used as a criterium of a 
low-temperature, equilibrated (Bose-condensed) system. 

(iii) Dynamical response of the non-equilibrium 
system (virith non-resonant pumping and photon 
losses) 

Below, we assume that relaxation mechanisms are not 
effective, and can not lead the polariton system to an 
equilibrium (thermal) state. The system is, thus, de- 
scribed by a density matrix, which is obtained from a 
master equation that takes care of photon losses through 
the cavity mirrors and incoherent (non-resonant) pump- 
ing. Details can be found in Ref. [2]. We solve the 
master equation in the stationary (t — > 00) limit in order 
to obtain the quasiequilibrium distribution, p^°°\ 

The response to the terahertz probe is computed in the 
linear approximation, where the probe does not modify 
the quasiequilibrium distribution. We adopt a compu- 
tational scheme similar to the one used for the photo- 
luminescence response [13]. The starting point is the 
first-order correlation function: 
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gdjj = (ii^ij - ^ij)gd,ij 
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with boundary conditions at i — > 00, t = 0: 



K 



(8) 



where, in the last step, we used the fact that p'^j is 
approximately diagonal in the energy representation (13j . 

In Eq. ([7|, K is the loss rate, 0.1 ps~^ in our model. 
The widths, F/j, are computed from: 
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where P is the pumping rate, and Nup{I) is the number 
of states with polariton number Npoi (/) + 1 used to solve 
the equations. 

The general solution of the linear system, Eq. ([t]), is 
written in terms of the eigenvalues, A„, and eigenvectors, 
of the matrix Bjj^mk defined by the r.h.s. of Eq. 



7L that is: 



gd.ijir) = ^C„exp(A„r)X}7, 
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written in terms of the auxiliary function: 



gd,ij^{{\J)mt + T)d{t)), 
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where | J) are A^porpolariton functions with total angular 
momentum L = 0, and the |/) are A'porpolariton func- 
tions with L = 1. Because of the Quantum Regression 
Theorem |15j . gd,ij satisfies the same equation as the 
density matrix, that is |14j : 



where the coefficients C„ are determined from the bound- 
ary conditions, Eq. 

The Fourier transform of Eq. ^ defines the response 
spectral function to the terahertz probe in the quasi- 
equilibrium system: 



I, J n 



where Dij_„ — ( J|d^|/)C„xjj\ and superscripts (r), (i) 
refer to the real and imaginary parts of the magnitudes, 
respectively. 

A simplified and more intuitive expression comes from 
the diagonal terms of Eq. ([7|. 14J Notice that, for ex- 
citation energies huj > 1 meV, the diagonal is at least 
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FIG. 3: (Color online) The lowest Npot = 2 states with L = 
and L = 1 in the model. Npoi — 1, L = states are drawn as 
a reference, and also in order to show that L = bands with 
different Npoi numbers are almost parallel. A big number of 
near zero-energy dipole transitions are possible in the Npoi = 
2 sector. 



10 times higher than the ofi-diagonal elements (because 
K — 0.1 ps"^). Neglecting the off-diagonal terms, we get: 



is 
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As compared with Sq, the non-equilibrium spectral 
function includes also contributions from the excited 
states, I J), which may have relatively high occupation 

probabilities, as can be seen, for example, in Fig. 

6 of Ref. [14]. On the other hand, the dipole matrix ele- 
ments for transitions originated in excited states could be 
much stronger than ground-state dipole elements. This 
statement follows from the energy-weighted sum rule for 
dipole transitions ^16j : 



J2^Eu\{I\d\J)f^C, 



(13) 



where constant C does not depend on the indices J. 

The sum in Eq. ( jlS] ) reduces to a single term when 
the oscillator strength from state | J) is concentrated on 
a single state, |/). Then, if there were excited states 
I J) for which the dominant transitions have AEjj ~ 0.1 
meV, for example, their contribution to Sne would be 
100 times stronger than the ground state contribution. 
This is, indeed, what one sees in the spectral function. 
Fig. [l] (b). An extra factor of around 20 comes from 
the number of excited states. We have drawn in this pic- 
ture the non-equilibrium spectral function for pumping 
rates, P, corresponding approximately to the same situa- 
tions depicted in Fig. [l](a). That is, the mean polariton 

number {{Npoi) = EjPjl^Npoi{J)) for P = 0.01 ps"!, 
for example, is around 4, etc. In Fig. |3j we show that 
near zero-energy dipole transitions are very common in 
our model, and should be very common also in micropil- 
lars with embeebed quantum wells because of the exciton 
near flat band. 

In conclusion, we expect the absorption spectral func- 
tion for a non-equilibrated polariton system to be peaked 
at near zero energies, in clear contrast with the Bose- 
condensed system, whose spectral function is peaked at 
the GDR. The dependence on Npoi is also very differ- 
ent. In the thermalized system absorption increases with 
increasing polariton number, whereas in the nonequilib- 
rium system it decreases as the pumping rate increases. 
Thus, terahertz absorption could be a sharp criterium al- 
lowing to discriminate between the thermalized and the 
non-equilibrium scenarios. 
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